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Abstract 

Utilizing recent DIS measurements ((Tr, -F2,3,l) and data on hadronic dilepton pro- 
duction we determine at NNLO (3-loop) of QCD the dynamical parton distributions 

of the nucleon generated radiatively from valencelike positive input distributions at 
an optimally chosen low resolution scale {Qq < 1 GeV^). These are compared with 
'standard' NNLO distributions generated from positive input distributions at some 
fixed and higher resolution scale (Qq > 1 GeV^). Although the NNLO corrections 
imply in both approaches an improved value of typically Xnnlo ~ 0.9xnlO' 
present DIS data are still not sufficiently accurate to distinguish between NLO re- 
sults and the minute NNLO effects of a few percent, despite of the fact that the 
dynamical NNLO uncertainties are somewhat smaller than the NLO ones and both 
arc, as expected, smaller than those of their 'standard' counterparts. The dynamical 
predictions for Fl{x, Q^) become perturbatively stable already at = 2 — 3 GeV^ 
where precision measurements could even delineate NNLO effects in the very small-x 
region. This is in contrast to the common 'standard' approach but NNLO/NLO dif- 
ferences are here less distinguishable due to the larger la uncertainty bands. Within 
the dynamical approach we obtain (M|) = 0.1124±0.0020, whereas the somewhat 
less constrained 'standard' fit gives as (Af|) = 0.1158 ±0.0035. 



1 Introduction 



Within the dynamical parton model approach the predicted small Bjorken-x behavior 
of structure functions is entirely due to QCD dynamics at x ^ 10^^. This is due to 
the fact that the parton distributions at Q"^ ^ 1 GeV^ are QCD radiatively generated 
from valencelike positive input distributions at an optimally determined low input scale 
Qo = /i^ < 1 GeV^ (where 'valencelike' refers to > for all input distributions 
/i^) oc — x)^f , i.e., not only the valence but also the sea and gluon input 

densities vanish at small x). Originally, its characteristic unique steep small-x predictions 
for the experimentally then unexplored region x < 10~^ [D 121 13] were subsequently first 
confirmed in |H |5]. With the advent of further high-precision data in recent years, the 
original dynamical parton distributions had to be updated [HI El IE] but the characteristic 
steep small-x behavior of the sea and the gluon distributions as x —>■ remained essentially 
very similar. 

Alternatively, in the common 'standard' approach, e.g. [9-21], the input scale is fixed 
at some arbitrarily chosen Qo > 1 GeV and the corresponding input distributions are less 
restricted. For example, the observed steep small-x behavior {aj < 0) of the gluon and sea 
distributions has to be fitted, allowing even for negative gluon distributions [IS HH [IH 
[TU] . i.e. negative cross sections like a negative longitudinal structure function Fl{x,Q'^). 
Furthermore the associated uncertainties encountered in the determination of the parton 
distributions turn out to be larger, particularly in the small-x region, than in the more 
restricted dynamical radiative approach where, moreover, the evolution distance (starting 
at Qo < 1 GeV^) is sizeably larger. 

In the present paper we extend our most recent LO and NLO dynamical analysis [S] 
to the next-to-next-to-leading order (NNLO) of QCD. For consistency reasons we only 
consider deep inelastic scattering (DIS) and Drell-Yan dimuon production data where all 
required theoretical NNLO ingredients are available by now, except the ones for heavy 
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quark production. High-pr hadron-hadron scattering processes will now not be considered 
since so far they are only known up to NLO. Furthermore we compare these 'dynamical' 
results for the radiativcly generated parton distributions arising from a valencelike pos- 
itive input at Qo < 1 GeV with the ones obtained from the common NNLO evolution 
approach being based on a 'standard' non-valencelike input at Qo > 1 GeV. In addition 
we shall analyze their associated uncertainties. Section 2 will be devoted to a discussion of 
some theoretical issues relevant for our NNLO analysis, in particular concerning our NNLO 
Q^-evolution algorithm for parton distributions, as well as of the relevant couplings and 
coefficient functions for relating them to the various neutral current structure functions 
required for calculating 'reduced' DIS cross sections. In Sect. 3 we present our quantita- 
tive results for structure functions, in particular our dynamical small-x predictions, and for 
hadronic dilepton production. The results related to the longitudinal structure Fl{x, Q^) 
are discussed in Sect. 4. Furthermore these dynamical results, together with their asso- 
ciated la uncertainties, are compared with the ones obtained from a common 'standard' 
approach. Our conclusions are summarized in Sect. 5. 



2 Evolutions of parton distributions and structure 
functions 



Our NNLO analyses will be performed within the modified minimal subtraction (MS) 
factorization and renormalization scheme. Heavy quarks (c, b, t) will not be considered as 
massless partons within the nucleon, i.e. the number of active (light) fiavors Uf appearing in 
the splitting functions and the corresponding Wilson coefficients will be fixed, Uf — 3. This 
defines the so-called 'fixed fiavor number scheme' (FFNS) which is fully predictive in the 
heavy quark sector where the heavy quark fiavors are produced entirely perturbatively from 
the initial light {u, d, s) quarks and gluons - in full agreement with present experiments. 
Furthermore, in the evaluation of the running strong coupling as{Q^) it is nevertheless 
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consistent and correct to utilize the standard variable Uf scheme for the /3-function [22]. 
Up to NNLO, as{Q'^) = as{Q'^)/4:7T evolves according to 

2 

daJdlnQ' = -J2f3ka':+' (1) 

k=0 

where Pq = 11 - 2nf/3, Pi = 102 - 38rij/3 and /?2 = 2857/2 - 5033n//18 + 325n^/54. 
Here we utilize the exact numerial (iterative) solution for aslQ"^) since such an accuracy 
is mandatory, in particular at NNLO, in the low region relevant for the valencelike 
approach [UIH]. Since Pk is not continuous for different nj, the continuity of as{Q^) requires 
to match ai"^^ at Q = rrih (h = c, b, t), i.e., a^J^^\mh) = a^J^^ ^\mh). These naive matching 
conditions get corrected at NNLO |23t 124^ 125] by a marg inalterm-(ll/727r2)[ai"^ ^\mh)f 
which will be neglected. Furthermore we have chosen rric = 1.3 GeV, = 4.2 GeV and 
rrit = 175 GeV, which turned out to be the optimal choices for our NLO analysis of heavy 
quark production [8j. 

The Mellin n-moments of the parton distributions f{x, Q^), 

f{n,Q')= [\xx'^-'f{x,Q') , (2) 







where f = q,q, g, evolve according to 



^^^Q2 = ^(««> ^m^^ Q ) (3) 

which refers to the coupled flavor-smg'/et evolution equation for q = (S, g)"^, S = "^giq+q), 
and 

P(a„ n) = Y: P^in) = E % % (4) 

fc=0 k=0 \ ^99 J 

with the well known LO and NLO splitting functions P/J''* and Pij \ respectively, and the 
NNLO (3-loop) P/j^^ have been calculated in p6]. Any obvious Q^- and/or n-dependence 
will be suppressed as far as possible. The 2x2 matrix evolution equation (3) can be 
formally solved recursively [271 128] with the result 

<f(n, Q2) = [L + asUiL - asoLUi + a^f/aL - a,a,of/i^f/i - a%L{U2 - f/i )] Ql) (5) 



where a^o = oisiQD/^'n- This is the so-called truncated solution where all redundant 0{a?g) 
terms are disregarded. The LO evolution operator L = L{as,aso,n), relevant for the LO 
solution gLo('^)Q^) = L{as,aso,n)q{n,Ql), can be written as 



L{as,aso,n) = { — ] = e_ ( — ) + e+ ( — ) (6) 



where Rq = Pq/Po and with the projection matrices e± being given by 

e± = T-^^ [Ro - X4] (7) 
where A_(A+) denote the smaller (larger) eigenvalue of Rq, 



^ [pS^ + n^? ± \/ip!S^-p^? + ^pSpp^ ] , (8) 



i.e., Rq = A_e_ + A+e+. Furthermore 

with -Rfc=i,2 = Rk+YnZl Rk-iUi and ^fc = Pk/Po-Yli=i PiRk-i/Po- We have not performed 
any required matrix multiplication in (5) analytically, since such a procedure did not 
reduce the required computer time of the subsequent numerical analysis. Rather, we have 
performed all required matrix multiplications entirely numerically, using the n-moments 

,(2) 



of the NNLO splitting functions [26] appearing in (4) together with the standard LO 



P,f and NLO P}p ones (see, e.g. [26]). The Bjorken-x space results q{x,Q'^) are finally 
obtained by performing numerically a contour integral around the singularities of q{n, Q"^) 
in (5) in the complex n-plane in the standard way (see, for example, [H [28l 129]). 

In the fiavoT-nonsinglet (NS) sector we have a simple (uncoupled) evolution equation 
which, in n-moment space, reads 

dqNs{n,Q^) , \ f n-2\ f,n\ 

^l^Q2 = ^Ns(as, n) gNs(^, Q ) (10) 
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where, similarly to (4), P^s{0's,n) = Y^I^q ci'^'^^ P^^ (n) which refers to the NS splitting 
functions P^g and P^s (s^^' example, [301 EI])- These splitting functions govern the 
evolution of the usual NS combinations of parton distributions with gNS referring to g^g^ 3 = 
— (i^, g^g g = + (i^ — 25=*", etc., where = g ± g, and gj^s = '^qiQ ~ q)- The NNLO 
splitting functions Pj^g'' have been given in [31] where the well known LO Pj^g* and NLO 
Pj^g ones can be found as well. Since no matrices are involved in the NS evolution equation 
(10), its solution can be easily inferred from the singlet solution (5) where now we have 
Ui = -Pi and 2U2 = -R2 - RiUi and thus 

q^is{n,Q^)= l-{as-a,o)Ri-^{al-alQ){R2-Rl)-aso{as-aso)Rl LqKs{n,Ql) (11) 

with Pfc = Pns//^o - YlLif^A-i/Po and L = L{as,aso,n) = (a^/a^o)"'^''- Again, 
QTSisix, Q^) is obtained from a numerical Mellin-inversion of gNs(^, Q^)- 

We have tested our singlet and nonsinglet evolution codes using the PEGASUS program 
[32] for generating the 'truncated' solutions together with the commonly used toy input of 
the Les Houches and HERA-LHC Workshops [331 El]- For 10"'^ < a; < 0.9 we achieved 
an agreement of up to four decimal places in most cases which is similar to the required 
high-accuracy benchmarks advocated in [33l IM] . 

As already mentioned at the beginning of this Section we employ for our analysis the 
FFNS and fix the number of active light flavors Uf = 3 in all splitting functions Pj^p 
and in the corresponding Wilson coefficients to be discussed below. In this factorization 
scheme only the light quarks {u, d, s) are genuine, i.e., massless partons within the nucleon, 
whereas the heavy ones (c, b, t) are not. This scheme is fully predictive in the heavy quark 
sector where the heavy quark flavors are produced entirely perturbatively from the initial 
light u, d, s quarks and gluons with the full heavy quark mass mc^h t dependence taken into 
account in the production cross sections - as required experimentally [351 ESI EZl EH], in 
particular in the threshold region. However, even for very large values of Q^, ^ ^cfe' 
these FFNS predictions up to NLO are in remarkable agreement [71 [8] with DIS data 
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and, moreover, are perturbatively stable despite the common belief that 'non-collinear' 
logarithms ln(Q^/m^) have to be resummed for h = c,b,t. 



This somewhat questionable resummation of heavy quark mass effects using massless 
evolution equations, starting at the unphysical 'thresholds' = m\, is persued in the so- 
called zero-mass 'variable flavor number scheme' (VFNS) where also the heavy quarks are 
taken to be massless partons within the nucleon with their distributions being generated, 
e.g. up to NLO, from the boundary conditions h{x, m\) = h{x, m^) = 0. Hence this factor- 
ization scheme is characterized by increasing the number Uf of massless partons by one unit 
at = m\ starting from n/ = 3 at = rnl, i.e., c{x,ml) = c(x,m^) = 0. The matching 
conditions are fixed by general continuity relations ^\ HO] at the respective 'thresholds' 

= m\. Thus the 'heavy' Uf > 3 quark distributions are perturbatively uniquely gen- 
erated from the Uf — 1 ones via the massless renormalization group Q^-evolutions (see, 
e.g. [10,15]; a comparative qualitative and quantitative discussion of the zero-mass VFNS 
and the FFNS has been recently presented in |ll]). Sometimes one uses an improvement on 
this, now known as the general-mass VFNS [IIl[I21[I31linill21llSlllllllSlllS], where mass- 
dependent corrections are maintained in the hard cross sections. This latter factorization 
scheme interpolates between the zero-mass VFNS and the (experimentally required) FFNS 
used for our analysis. 

In order to avoid any further dependence on model assumptions, we choose to work with 
experimentally directly measurable quantities, as has been done in P], like the 'reduced' 
DIS one-photon exchange cross section a,. = F2 — {y"^ /Y+)Fl together with the full neutral 
current (NC) cross sections |17] 




-1 




dx dQ"^ 
^ and 



2 



.NC 



(12) 



where a = 1/137.036, Y± = 1 ± {1 - y) 




F,^,-v^KFll + {vl + aiyFi^, 
—agnF^^ + 2veaen'^F^ , 



(13) 
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with Ve = —I + 2sin^^i4/, Oe = — I and = Asin^ 6w cos^ OwiQ"^ + Ml)/Q'^, using 
sin^^vi/ = 0.2312 and Mz = 91.1876 GeV. As in our previous NLO analysis [S] it turned 
out, however, that fitting just to the usual (one-photon exchange) F2{x,Q^) gives rather 
similar results. Defining J-'2,l = ^^2,1/^, the n- moments in (2) of these structure functions 
^2,l{x,Q'^) and F^^XjQ"^) can, for nj = 3 light fiavors, be written as 

F^^ = C'3,Ns(a"gNS + «3 W + (14) 
where af = - a^), 4 = + 4 " 2a^) and a± = | Eg=„,d,s with 

a+ = el-2eqVeVgK+ {vl + aDK^ 

a~ = —2eqaqK + AveaeVgagK^ (15) 

where = ±^ — 2egSin^6'i^ and = ±^ with ± according to whether q is a u- or 
d-type quark. The Wilson coefficients are generically expanded as [30l C2,3{as,n) = 
Sfc=o '^sC2^3(^) where at LO C2°ns = ^2°] = cf^^ = 1 and C2°] = 0. The appropriate 
NLO (1-loop) coefficients be found, for example, in |3U]. The NNLO (2-loop) 

coefficients C2 3 have been originally calculated in [501 EI] and, for definiteness, we take 
c^2Ns = '^2^NS' '^2^9 ^^"^ '^2^] from [S2], and Cg^^g = Cg^^g from [51j. Since the longitudinal 
structure function Fl = F2 — 2xFi vanishes at LO, it has become common [50] to consider 
the first nonvanishing 0{as) contribution to Fl as the LO one, i.e., the perturbative 
expansion up to NNLO now reads CL{as,n) = Yll:=i ^s'^^l\''^) ■ "^^^ relevant NNLO (3- 
loop) coefficients c^^ have been calculated in [521 153], where also the well known LO c^^^ 
and NLO c)^ coefficients can be found. Although we perform all calculations in Mellin 
n-moment space, it should be nevertheless mentioned that in Bjorken-x space the simple 
products in (14) turn into the standard convolutions of the Wilson coefficients with the 
parton distributions. 

In the medium to large x-region the relevant kinematic nucleon target mass (TM) 
corrections are also taken into account for the dominant 'light' F2 structure function in (14) 
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(with 'light' referring to the common u, d, s (anti-)quarks and gluon initiated contributions) 
according to 



F2,TM(ri,Q') = I rfxx"-2F2,TM(x,Q') 



£=0 



\ Q2 J ^\{n - 2)! (n + 2l){n + 21-1) \\Q^ ) I ^ ' 



where higher powers than {m\/Q'^Y ^"^^ neghgible for the relevant x < 0.8 region, as can 
straightforwardly be shown by comparing (16) with the well-known exact expression in 
Bjorken-x space 



So far we have discussed only the contributions of light partons (u, d, s, g) to structure 
functions in (14), -Fj'=2L3- The total structure functions Fi{x,Q'^) = Fj^^^^ _|_^heavy j,gq^j],g 
also the knowledge of the (subleading) heavy quark contribution F^'^^^^y = + Fj' at 
fixed-order of perturbation theory. (Top quark contributions are negligible.) The LO 
0{as) contributions to F2L, due to the subprocess 'y*g — > hh with h = c,b, have been 
summarized in [6j, and the NLO 0{a1) ones are given in [551 EE]. The NNLO 0{af) 
3-loop corrections to F^ and first rudimentary contributions to F2 have been calculated 
recently EHl EH] for ^ ml, but these asymptotic results are neither applicable for 
our present investigation nor relevant for the majority of presently available data at lower 
values of Q^. Our ignorance of the full NNLO 0{a^) corrections to constitute the 
major problem for any NNLO DIS analysis in the FFNS. (It should be mentioned that 
we have attempted to mimick the NNLO contributions by naively assuming them to be 
down by one power of times the NLO terms multiplied by a constant i^"-factor, but the 
fit results were insensitive to such an ad hoc correction. However, this approach (guess) 
appears to be not appropriate since playing the same game at NLO, i.e., times LO times 
a K-factor can not reproduce the correct NLO results in the relevant kinematic region of 
X and (5^.)n Therefore, the heavy flavor contributions F2X, are taken as given by fixed 



further (inconsistent) 'check' of the relevance of the unknown massive NNLO coefficient functions 
can be made by comparing the predicted charm and bottom structure functions using our new NNLO 
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order NLO perturbation theory [551 156] as in our previous more restricted NNLO analyses 
[nni EI] • This is also common in the literature fTEl \W\ [T7] and the error in the resulting 
parton distributions due to NNLO corrections to heavy quark production is expected [15] 
to be less than their experimental errors. These contributions are gluon g{x, dominated 
and the factorization scale, also of the remaining parton distributions, should preferably 
chosen [62] to be iJ?p = 4:m\, although a much larger choice like = 4(Q^ + 4m^) leaves 
the NLO results essentially unchanged [TJ [H|. The NNLO heavy quark contributions to 
F^'-^ are not known either, but here they vanish in LO and are already negligibly small in 
NLO at the relevant large values of as discussed in [8]. 

More recently the NNLO corrections to the rapidity distribution (Pa/dM dy of Drell- 
Yan (DY) dilepton production of mass M has been calculated as well [631 El] • This allows 
to include DY data as well for performing a fully consistent analysis up to NNLO. Needless 
to say that the DY pp and pd dilepton production data are instrumental in fixing d — u (or 
d/u). Only the usual high-p^ inclusive jet production data of hadron-hadron scattering 
have to be disregarded where the NNLO corrections have not yet been calculated. The 
LO and NLO corrections to the DY process are well known (for a summary, see, e.g., [65] ) 
and for our full NNLO analysis we used the routine developed in [66| based on the results 
of [631 El- 

Finally, the evaluation of the uncertainties of our NNLO parton distributions is per- 
formed in the same way as of our recent NLO ones [H] which followed the line of [HZl EHl EH]- 
The uncertainties Aoj = — a° of the central free fit parameters a°, corresponding to the 
minimal Xo^ constrained by A^^ < with the tolerance parameter T chosen to be 

= = v^/(1.65)2 - (4.5)^ i.e., T being slightly smaller than in [8] due to the 

smaller total number of data points considered, = 1568, because the high-py jet data 

(anti)quark and gluon distributions, as presented in Sec. 3, and the currently known massive NLO coef- 
ficient functions with the fully consistent NLO predictions (e.g. jS]) based on NLO parton distributions. 
These predictions turn out to be indistinguishable, except at very small values of x, x < 10~^, where the 
'NNLO results' are about 10-15% smaller than the NLO ones as shown for example in Figs. 8 and 9 of |8]. 
This is still fully consistent with the cc and bb HERA DIS data at very small x [5511551 [571 [55] . 
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and the DIS data for semi-inclusive cc- and 66-production cannot consistently be included 
in a global NNLO fit for the time being. 

3 Quantitative results and dynamical small-x predic- 
tions 

Now we extend our recent dynamical LO and NLO(MS) analysis [S] to NNLO (MS). The 
valencelike input distributions xf{x, Ql) at the input scale Qo = fi < 1 GeV, referring to 
the flavor nonsinglet (valence) densities m„, d^, A = d — u and to the valencelike densities 
d + u, s = s and g in the singlet sector, are generically parametrized as 



Since the data sets we are using are insensitive to the specific choice of the strange 
quark distributions, we continue to generate the strange densities entirely radiatively 
[Zl [8] starting from s{x, Q^) = s{x, Q^) = in the dynamical valencelike approach 
where Qq < 1 GeV. For comparison we also study the common standard evolution ap- 
proach, being based on a non-valencelike input at Qo > 1 GeV, where we choose as usual 
s{x, Qg) = s{x, Ql) = [u{x, Ql) + d{x, (5o)]/4. Furthermore, since all our fits did not require 
the additional polynomial in (17) for the gluon distribution, we have set Ag = Bg = 0. This 
left us with a total of 21 independent fit parameters, including as- As suggested in and 
done in our previous NLO analysis [8J, we included in our final error analysis only those 
parameters that are actually sensitive to the input data set chosen, i.e. those parameters 
that are not close to 'flat' directions in the overall parameter space. With current data, and 
our functional form (17), 13 such parameters, including a^, are included in our final error 



xfix,Ql) = Nfx"^ {1 - xp {1 + Af^Jx + Bfx) , 




subject to the constraints u^dx = 2, d^dx = 1 and 
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analysis. The remaining highly correlated ill-determined eight polynomial parameters Af 
and Bf, with uncertainties of more than 50%, were held fixed. 

These free parameters have been fixed using the following data sets: the HERA ep 
measurements jTOl [TH [72l [731 El] for > 2 GeV^ for the 'reduced' cross sections 
and ar,NC in (12); the fixed target F| data of SLAG [73], BCDMS ^76], E665 [77] and 
NMC [78], subject to the standard cuts > A GeV^ and = ^^(1 _ i) + ^2 > 
GeV^ together with the structure function ratios F^/F^ of BCDMS [79], E665 [HQ] and 
NMC [81]. Furthermore the Drell-Yan muon pair production data of E866/NuSea [82] for 
d'^aP^/dM dxF with N = p,d have been used as well as their asymmetry measurements [HS] 
for a^'^/a^P. The DY data are always given in terms of Xi?-distributions, whereas the NNLO 
expressions have been given in terms of the dilepton rapidity ^-distributions [631 EH 166] . 
Since experimentally the dilepton px is small (below about 1.5 GeV) as compared to the 
dilepton invariant mass M ^ 5 GeV, we have checked that it can be safely neglected and 
the two distributions can be related using leading order kinematics, as has been done in 
[T7] : d'^a/dMdxp = {xi + XiY^d^a / dM dy where Xi,2 = (MVs)^/^e=^^, Xp = Xi - X2 and 
xi+ X2 = (xjp + AM"^ / sy/'^. All these data sets correspond to 1568 data points. 

The parameters obtained from our NNLO dynamical fit for the input distributions at 
the optimal input scale Ql = /Unnlo ~ O-^^ GeV^ are given in Table 1, and the ones 
for our standard fit, corresponding to the choice Ql = 2 GeV^, in Table 2. Due to the 
fact that the higher the perturbative order the faster ^^(Q^) increases as decreases, 
a NNLO analysis is expected to result in a smaller value for asiM"^) than a NLO fit in 
order to compensate for this increase. Our dynamical analysis indeed yields a smaller 
value at NNLO in Table 1, a^^L°(M|) = 0.1124 ± 0.0020, as compared to our NLO fit 
[8] which resulted in Q;f^°(M|) = 0.1145 ± 0.0018. Both values lie, however, within a 
la uncertainty. (Similar results were obtained in a previous dynamical fit [HI] which was 
performed for a restricted set of (mainly small-x) DIS data.) The same holds for our 
NNLO 'standard' fit result in Table 2, a^^^^{M^) = 0.1158±0.0035, to be compared with 
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the NLO result [8] of 0.1178 ± 0.0021. Our standard NNLO fit result for as(M|) in Table 
2 agrees, within errors, with the results of [13 US] and [IT] , 0.1143 ± 0.0014 and 0.1128 ± 
0.0015, respectively, and is compatible with the one obtained from a 'standard' fit [60j to 
a restricted set of small-x DIS data. Data for high-p^ jet production in hadron-hadron 
scattering should not be included in a consistent NNLO analysis, since such processes are 
theoretically known only up to NLO. Including them nevertheless in a standard NNLO 
analysis requires [HI El] generally larger values of as(M|). Previous standard NNLO fits 
{Qo > 1 GeV) considering only the fiavor non-singlet (NS) valence sector of structure 
functions [HHl IHEl EZ] resulted in somewhat smaller values of as(M|) than in Table 2 but 
remain within a Icr — 2(T uncertainty. A similar NS valence analysis [88j as well as a full 
analysis [H9] being based, however, on incomplete calculations of the moments of 3-loop 
anomalous dimensions (splitting functions) yielded slightly larger values of at NNLO, 
as(M|) ~ 0.117, with estimated errors large enough so as to comply with our result in 
Table 2. For a more detailed and comparative recent discussion of NLO and NNLO results 
the interested reader is referred to [9^. In general, the NNLO fits result in a better 
(smaller) than the NLO ones 0, typically Xnnlo — 0-9Xnlo- 

It should be noticed that our a^-uncertainty in Table 2 is about twice as large as the one 
obtained in a comparable standard NNLO analysis [17] where the high-p^ jet data have 
been disregarded for consistency reasons as well. Without these data the gluon distribution 
is little constrained in the medium to large x-region where it plays an important role for 
the Q^-evolution at small values of x due to the convolution with the dominant Pgg\ This 
a^-uncertainty remains sizeable irrespective of the choice of the input scale Qq > 1 GeV. 
Only within a Bayesian treatment of systematic errors, by taking into account point-to- 
point correlations [HH IS21 [S3], the uncertainty of as(M|) turns out to be about two times 
smaller [T5l [T6l [T71 [9T| [92| [93] . On the other hand the a^-uncertainty of our dynamical 
fit in Table 1 is also about half as large as the 'standard' one in Table 2. Apart from the 
larger evolution distance, this is due to the strongly constrained valencelike input gluon 
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distribution xg{x,Ql) = NgX°'^{l — at = 0.55 GeV^ in the small-x region where 
ttg ~ 1, according to Table 1; consequently the energy-momentum sum rule (18) sufficiently 
constrains xg{x, Qq) in the medium to large x-region as we shall see below. 

Our dynamical NNLO valence and valencelike (sea and gluon) input distributions at 
Qo ~ A^NNLO ~ 0-^^ GeV^ are shown in Fig. 1, according to the parameters in Table 1, 
together with their la uncertainties. A comparison with our previous NLO results [8] 
(dashed curves) shows that and are now somewhat enhanced around x = 0.1 to 0.2 
and that a strong and clear valencelike small-x behavior of the NNLO gluon input is now 
required, ag = 0.994±0.379, as compared to ag ~ 0.5168± 0.4017 at NLO [Hj. Furthermore 
there is also a strong enhancement of the NNLO gluon over the NLO one around x = 0.1 
and a sizeable depletion at larger values of x. The valence distributions of our standard 
analyses {Qq = 2 GeV^) are compared with the standard NNLO ones of Alekhin, Melnikov 
and Petriello [T7j (AMP06, Ql = 9 GeV^) and the standard pure NS analysis of Bliimlein, 
Bottcher and Guffanti [87J (BBG06, = 4 GeV^) in Fig. 2 at = 4 GeV^. In the 
relevant valence x-region, x ^ 0.1, we conffim the NNLO BBG06 results, in particular the 
enhancement of xd^ with respect to the NNLO result of AMP06. In any case we, as well 
as BBG06, observe a significant enhancement of the NNLO xu^ and xd^ with respect to 
the NLO results. The valence distributions are very robust with respect to the choice of 
the input scale Ql since our dynamical valence distributions at = 4 GeV^ practically 
coincide with the standard ones shown in Fig. 2. 

The distinctive valencelike gluon input at low Q"^ < 1 GeV^ in Fig. 1 implies a far 
stronger constrained gluon distribution at larger values of as compared to a gluon 
density obtained from a 'standard' fit with a conventional non-valencelike input at > 1 
GeV^, Qq = 2 GeV^ say, as can be seen in Fig. 3. In contrast to the standard NLO 
results [H] (shown by the dotted curves in Fig. 3) the standard NNLO gluon input at 
Qo~Q^~'^ GeV^ is very weakly constrained at small x {ag = 0.0637 ± 0.1333, cf. Table 
2) and therefore (18) cannot sufficiently constrain it at larger values of x since, moreover. 
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high-pT jet data have not been taken into account for consistency reasons. Notice that this 
common standard NNLO input gluon distribution at Qq = 2 GeV^ is also compatible with 
a valencelike small-x behavior (a^ > 0) — a tendency already observed in |13j — and that 
our dynamical NNLO gluon distribution in Fig. 3 (solid curves) remains valencelike even 
at = 2 GeV^ (i.e. decreases with decreasing x). This is mainly caused by the NNLO 
splitting function Pgg in (4) which is negative and more singular in the small-x region [26] 
than the LO and NLO ones: for Uf = 3, xP^g\x) ~ -3147.66 In ^ + 14737.89 as x ^ 0, 
whereas xPgg\x) ~ 12 and xPgl\x) ~ —81.33. The uncertainties generally decrease as 
increases due to the QCD Q^-evolutions [67l|69j, but the ones of the dynamical predictions 
in the small-x region remain substantially smaller than the uncertainties of the common 
'standard' results as exemplified in Fig. 3. Furthermore, it is a general feature of any 
NNLO gluon distribution in the small-x region that it falls below the NLO one as can be 
seen in Fig. 3 by comparing the solid curves with the long-dashed ones, and the dashed- 
dotted curves with the dotted ones. In the first dynamical case the NNLO predictions for 
X < 10~^ are several-cr below the NLO ones. For comparison we also display in Fig. 3 the 
standard NNLO results of AMP06 [H] based on an input scale Ql = 9 GeV^. 

The dynamical sea distribution x{u + d) derives from a less pronounced (a^^+j < ag) 
valencelike input in Fig. 1 which vanishes very slowly as a; — > (a^_,_j = 0.1374 ± 0.0501, 
cf. Table 1). This implies that the valencelike sea input is similarly increasing with de- 
creasing X down to X ~ 0.01 as the sea input obtained by the common standard fit where 
a^i+j = —0.1098 ± 0.0122 according to Table 2. Therefore, the la uncertainty bands of our 
dynamically predicted sea distributions at larger values of in Fig. 4 are only marginally 
smaller than the corresponding ones of the standard fit. In contrast to the evolution of 
the gluon distribution in Fig. 3, the NNLO sea distributions in Fig. 4 lie always above the 
NLO ones in the small-x region, x ^ 10~^, and at not too large values of Q^. Here all 
NNLO sea distributions are rather similar, including the 'standard' one of AMP06 [IT] . 

A representative comparison of our dynamical and standard NNLO results with the 
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relevant HERA(H1,ZEUS) data on the proton structure function is presented 

in Figs. 5 and 6. It should be reemphasized that due to our valencelike input, the dy- 
namical small-x results {x ^ 10"^) are predictions being entirely generated by the QCD 
Q^-evolutions. This is in contrast to a common 'standard' fit where the gluon and sea 
input distributions in (17) do not vanish as x — > 0. For comparison we also display our 
dynamical NLO results [8] shown by the dashed curves. In all cases the data in Figs. 5 and 
6 are well described throughout the whole medium- to small-x region for ^ 2 GeV^ and 
thus perturbative QCD is here fully operative. At < 2 GeV^ the theoretical results fall 
below the data in the very small-x region; despite the fact that the NNLO results in Fig. 5 
are closer to the data than the NLO ones, this is not unexpected for perturbative leading 
twist-2 results since nonperturbative (higher twist) contribution^ to F2{x, Q^) will eventu- 
ally become relevant, even dominant, for decreasing values of Q^. Although the inclusion of 
NNLO corrections imply an improved value of x^, typically Xnnlo — O-O^Xnlo according 
to [H] and Tables 1 and 2, present high precision DIS data are not sufficiently accurate to 
distinguish between the NLO results and the minute NNLO effects of a few percent. This 
is illustrated in Fig. 7 where the experimental (statistical and systematic) errors are far 
bigger than the differences between the NLO and NNLO results. It should, however, be 
noticed that the NNLO Icr uncertainty band is somewhat narrower (reduced) than the one 
at NLO. The results are similar for our 'standard' fits. It has already been noticed that, by 
analyzing only the fiavor non-singlet valence sector of structure functions, NNLO effects 
cannot be delineated by present data in the medium- to large-x region, and moreover, un- 
certainties of NLO and LO analyses (such as higher twists, different factorization schemes 
and QED contributions to the QCD Q^-evolutions) turn out to be comparable in size to 
the NNLO 3-loop contributions [SB] . 

As already pointed out, the measurements of Drell-Yan dilepton production in pp and 

^Since operators of different twists do not mix under renormalization group evolutions, possible higher- 
twist effects do not influence the determination of the input parameters at the low input scale Qq — 
0.55 GeV^ relevant for the dynamical (valencelike) leading twist~2 distributions, with the latter being 
determined from data at > 2 GeV^ and > 4 GeV^ with W'^ > 10 GeV^ as discussed above. 
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pd collisions [821 [83] are instrumental in fixing A = d — u (or d/u) [M]. In Figs. 8 
and 9 we display our dynamical NNLO and NLO results, together with the ilex un- 
certainties, for the differential dimuon mass distributions for various average values of 
xp = Xi — X2 for pp and pd collisions, respectively. The 'standard' fit results differ only 
marginally. In the relevant kinematic region where high-statistics data exist, all three 
NNLO and NLO results shown agree within la. In Fig. 10 we show the result for the 
ratio aP'^/2aPP relevant for the DY asymmetry Ady = (c^^^ - aP'')/{(rPP + (jP"). Notice 
that (T^^ = d^a^^ /dMdxp oc d s [^(^1)^(^2) + (l{,X2)q{xi)\ in LO at a scale Q"^ = M^, 
where xi and X2 refer to the fractional momenta of the quarks in the beam (p) and the 
nucleon target {N), respectively. Experimentally xp > ^ {xi > X2) and consequently the 
Drell-Yan cross section is dominated by the annihilation of a beam quark with a target 
antiquark. 

4 The longitudinal structure function Fl{x^ Q'^) 

As discussed in Sect. 2 we have explicitly used for our analysis the experimentally directly 
measured 'reduced' DIS cross sections (12) which, for not too large values of Q^, are 
dominated by the one-photon exchange cross section = -F2 — {y'^ /Y^)Fl where y = 
Q'^/xs. The importance of using this quantity has been emphasized in [95]: the effect of 
Fl becomes increasingly relevant as x decreases at a given Q^, where y increases. This is 
seen in the data as a flattening of the growth of ar{x, Q"^) as x decreases to very small values, 
at fixed Q^, leading eventually to a turnover (cf. Fig. 11). At lower values of in Fig. 11 
ist was not possible in [95] to reproduce this turnover at NLO. This was mainly due to 
the negative longitudinal cross section (negative Fl(x, Q^)) encountered in [131 195]. Since 
all of our cross sections and subsequently all structure functions are manifestly positive 
throughout the whole kinematic region considered, our dynamical NLO [8] and NNLO 
results in Fig. 11 are in good agreement with all small-x HERA measurements [70-74]. 
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The same holds true for our 'standard' NLO [8] and NNLO results which, besides having 
slightly wider uncertainty bands, are almost indistinguishable from the dynamical ones 
shown in Fig. 11. In Fig. 12 we display our NNLO results for at different proton 
beam energies i?p, relevant for most recent HI measurements [96], where the turnover at 
small X becomes more pronounced at smaller energies because of the larger values of y. 
Our dynamical small-x predictions are fully compatible with the (preliminary) HI data 
presented in [91) j . 

Turning now to itself we note that the ra-moment equation (14) for becomes in 
Bjorken-x space 

x-^F^{x, Q') = Cl,ns ® g^5,3 + ^ <lNS,s) + I (Cl,, ® S + C^,, ®g)+ x~'F^''''^x, Q') 

(19) 

where (8> in the light parton sector denotes the common convolution, and the weak 
contributions in (13) and (15) have been neglected for <^ M| relevant for our present 
interest. We have also added the heavy quark (charm, bottom) contribution F^^^"^^ for 
which we use, as discussed in Sect. 2, the NLO(2-loop) expressions also in NNLO due to our 
ignorance of the 0{al) NNLO heavy quark corrections. (This inconsistency is here of minor 
importance since F£ - and even more so - is a genuinely subdominant NLO contribution 
to the total Fl, which holds of course also at LO (cf. Figs. 13 and 14)). Following the 
notation of Sect. 2, the perturbative expansion up to NNLO of the coefficient functions in 
(19) reads Cz.,(«.,x) = ELi In LO, = fx, c^, = 0, c^^l = 2Ax{l - x) 

and the flavor-singlet quark coeefficient function is decomposed into the non-singlet and 
'pure singlet' contribution, c^^^ = c^l'ms^ ^^I'ps- Sufficiently accurate simplified expressions 
for the NLO [971 [981 |99] and NNLO [Sg coefficient functions c^l and cfl, respectively, 
have been given in [53]. It has been furthermore noted in that especially for Ci^g 
both NLO and NNLO contributions are rather large over almost the entire x-range. Most 
striking, however, is the behavior of both singlet coefficient functions g and C^^g in 
(19) at very small values of x\ the vanishingly small LO parts {xc^l\ ~ x^) are negligible 
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as compared to the negative constant NLO 2-loop terms, which in turn are completely 
overwhelmed by the positive NNLO 3-loop singular contribution xc"^\ ~ In |. This latter 
singular correction might be indicative for a perturbative instability at NNLO [53j but it 
should be kept in mind that a small-x information alone is insufficient for reliable estimates 
of the convolutions occuring in (19) when evaluating physical observables. 

Our dynamical LO, NLO and NNLO predictions for the total Fl are displayed in Fig. 13, 
together with the small subdominant charm contributions at LO and NLO. These predic- 
tions become perturbatively stable already at ^ 3 GeV^ where the gluon contribution 
becomes dominant and where precision measurements could even delineate NNLO effects 
in the very small-x region. It should be noted that the error bands at smaller values of 
(where the gluon and quark contributions are still comparable) are rather small. This is 
caused by compensating effects in the error analysis, although the individual light quark 
and gluon contributions have wider error bands. Within the common 'standard' approach 
the absolute values of the NNLO and NLO results in Fig. 14 differ by more at smaller 

than the dynamical ones, but the differences between the NNLO and NLO results are 
here less distinguishable due to the larger la uncertainty bands which partly overlap in 
the very small-x region. It should, however, be noted that it is somewhat deceptive to 
compare the error bands in Figs. 13 and 14: in the 'dynamical' approach (Fig. 13) the 
errors are strongly reduced due to the evolution from about 0.5 GeV^, in contrast to the 
much smaller evolution distances in the 'standard' approach (Fig. 14) where the evolution 
starts at the input scale = = 2 GeV^. It should furthermore be noticed that the 
NLO/NNLO instabilities implied by the standard fit results obtained in [T31 ES] at 
^ 5 GeV^ are far more violent than the ones shown in Fig. 14 which is mainly due to 
the negative longitudinal cross section (negative Fl{x,Q'^)) encountered in [131 ES]- The 
perturbative stability in any scenario becomes in general better the larger Q^, typically 
beyond 4 — 5 GeV^ [131 EBl EU |95], as evident from Figs. 13 and 14. This is due to the 
fact that the Q^-evolutions eventually force any parton distribution to become sufficiently 
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steep in x. It should be mentioned that the sizeable discrepancies between NNLO and NLO 
predictions at = 2 GeV^ and x ~ 10^^ in Figs. 13 and 14 are not too surprising since 
~ 2 GeV^ represents somehow a borderline value for the leading twist-2 contribution 
to become dominant at small-x values. This is further corroborated by the observation 
that the dynamical NNLO and NLO twist-2 fits slightly undershoot the HERA data for 
F2 at 2 GeV^ in the small-x region (cf. Fig. 5), which indicates that nonperturbative 

(higher twist) contributions^ to structure functions become relevant for Q"^ ^ 2 GeV^ 13 E]- 

For completeness we finally compare in Fig. 15 our NNLO dynamical and standard 
(leading twist) predictions for Fl(x, Q^), together with their ±lcr error bands, with a 
representative selection of (partly preliminary) HI data [721 E31 HU^ llOlj llU2j at fixed 
W ~ 276 GeV. For comparison we also show in Fig. 15 our NLO results [8] which have la 
uncertainty bands similar to the NNLO ones. All our NNLO and NLO results for F^, being 
gluon dominated in the small-x region, are in full agreement with present measurements 
which is in contrast to expectations [12], [13], [95] based on negative parton distributions and 
structure functions at small values of x. To illustrate the manifest positive definiteness of 
our dynamically generated structure functions at > /i^ (/Wnnlo — 0-^5 GeV^, /iNLo — O-^ 
GeV^), we show Fl{x, Q"^) in Fig. 15 down to small values of although leading twist-2 
predictions need not necessarily be confronted with data below, say, 2 GeV^. 



5 Summary and conclusions 



Utilizing recent DIS structure function measurements (-F2,3,l and the 'reduced' cross section 
(Tr) and hadronic Drell-Yan dilepton production data, our previous LO and NLO global 
fit analyses for the dynamical parton distributions of the nucleon [BJ have been extended 
to NNLO of perturbative QCD. The small-x {x ^ 10~^) structure of dynamical parton 
distributions is generated entirely radiatively from valencelike, manifestly positive, input 
distributions at an optimally chosen input scale Ql < 1 GeV^. The NNLO predictions are 
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perturbatively stable with respect to the NLO ones and are in agreement with all present 
measurements for ^ 2 GeV^. In general, the NNLO corrections imply an improved value 
of typically Xnnlo — ^-^ Xnlo- Having augmented our analyses with an appropriate 
uncertainty analysis, it turned out that present DIS precision data are still not sufficiently 
accurate to distinguish between NLO results and the minute NNLO effects of a few percent, 
despite of the fact that the dynamical NNLO la uncertainties are somewhat smaller than 
the NLO ones. Inclusive high-p^ jet data were disregarded for consistency reasons, since 
NNLO corrections have not yet been calculated. Nevertheless the valencelike input gluon 
distribution remains sufficiently constrained, via the energy-momentum sum rule, also 
at larger values of x (which is in contrast to a common 'standard' fit approach). It is 
interesting to note that our dynamical NNLO gluon distribution remains valencelike even 
at (5^ = 2 — 3 GeV^ (i.e. decreases with decreasing x, cf. Fig. 3) which is mainly caused 

(2) 

by the dominant NNLO gluon-gluon splitting function Pgg' which is negative and more 
singular than the LO and NLO ones, Pgg {x) ~ — ^ In ^. The drawback of any 

precision NNLO analysis at present is that the experimentally required heavy quark mass 
effects of heavy quark (charm, bottom) contributions can only be taken into account up 
to NLO because of our ignorance of the full NNLO C(q;^) corrections. 

Our dynamical distributions and predictions have also been compared with conven- 
tional ('standard') NNLO ones obtained from non-valencelike positive definite input dis- 
tributions at some arbitrarily chosen higher input scale Ql > 1 GeV^. For this purpose 
we have performed a common 'standard' fit as well, assuming Q"^ = 2 GeV^ (notice that, 
contrary to the dynamical approach, the finite small-x behavior of the input gluon and 
sea distributions is here fitted, and not dynamically generated by QCD evolutions). As in 
the dynamical approach, the NNLO corrections imply here also an improved x^, typically 
Xnnlo — 0-9Xnlo- '^^^ uncertainties of these less constrained 'standard' distributions 
are, as expected, larger than those of their dynamical counterparts, in particular in the 
small-a; region. 
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Our predictions for the longitudinal structure Fl^x^Q"^) and results for the 'reduced' 
DIS cross section (Tr(x, Q^) are in agreement with all HERA data and most recent HERA- 
HI measurements, in particular in the small-x region and down to = 2 GeV^. The 
dynamical NNLO/NLO predictions for Fl{x,Q'^) become perturbatively stable already at 

= 2 — 3 GeV^ where future precision measurements could even delineate NNLO effects 
in the very small-x region (the NNLO la uncertainty bands are here smaller than the NLO 
ones). This is in contrast to the common 'standard' approach but NNLO/NLO instabilities 
and differences are here less distinguishable due to the much larger la error bands. 

The strong coupling obtained from our dynamical NNLO analysis is af'^^*~'(M|) = 
0.1124 ± 0.0020 to be compared with a^^^{Ml) = 0.1145 ± 0.0018 at NLO [8J. The less 
constrained 'standard' approach at NNLO resuhed in a^^^°(M|) = 0.1158 ± 0.0035 to be 
compared with a^L°(M|) = 0.1178 ± 0.0021 at NLO 0. 

For our analysis in the 'fixed flavor number scheme' with Uf = 3 active light {u,d,s) 
flavors, we have developed our own (entirely numerical) NNLO Q^-evolution algorithm 
which has been tested by reproducing the appropriate Les Houches and HERA-LHC high- 
accuracy benchmarks |33l[3l]. A FORTRAN code (grid) containing our NNLO dynamical 
and 'standard' light (m, d, s; g) parton distributions can be obtained on request or directly 



from http : / /doom . physik . uni-dortmund . de/ pdf server 
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NSLO (MS) 




Uv 


dy 


d — u 


u + d 


9 


N 


4.4049 


13.824 


8.6620 


1.2316 


23.034 


a 


0.7875 


1.1778 


1.2963 


0.1374 


0.9940 


b 


3.6857 


5.6754 


19.057 


10.843 


6.7892 


A 


-1.1483 


-2.2415 


-6.8745 


-4.5634 




B 


4.5921 


3.5917 


19.402 


11.940 




xVdof 


0.986 (0.904) 


a.(M|) 


0.1124 ± 0.0020 



Table 1: Parameters of our dynamical input distributions as parametrized in (17) referring 
to an input scale of Ql = /^nnlo = 0.55 GeV^. Since the gluon distribution turned out to be 
insensitive to the polynomial terms in (17), we have set them to zero {Ag = Bg = 0). The 
total numbers of degrees of freedom is dof = 1568—21 = 1547. The x^/dof in brackets refers 
just to the DIS data where dof = 1178 - 21 = 1157. Furthermore ^^(//^nlo)/^ = 0.1522. 





NSLO (MS) 




Uv 


dx, 


d — u 


u + d 


9 


N 


3.2350 


13.058 


8.1558 


0.4250 


3.0076 


a 


0.6710 


1.0701 


1.1328 


-0.1098 


0.0637 


b 


3.9293 


6.2177 


21.043 


10.341 


5.4473 


A 


-0.5302 


-2.5830 


-7.6334 


-3.0946 




B 


3.9029 


3.8965 


20.054 


11.613 




xVdof 


0.947 (0.873) 




0.1158 ± 0.0035 



Table 2: As Table 1 but for the input parameters in (17) of the standard fit at an input 
scale Ql = 2 GeV\ where as{Ql)/TT = 0.1072. 
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X 



Figure 1: The dynamical NNLO valence {uv,dy,d — u) and valencelike {g^u + d) input 
distributions together with their ilcr uncertainties at = /^nnlo ~ ^-^^ GeV^. The 
central curves follow from (17) with the parameters given in Table 1. The strange sea 
s = s vanishes at the input scale. Our dynamical NLO input [8J at Ql = /Inlo — 0-^ GeV^ 
is also shown by the dashed curves for comparison. The la uncertainties at NLO [8] are 
comparable to the ones shown at NNLO, except for the NLO gluon at x > 0.3 which is 
stronger constrained due to the light high-p^- jet data [8j. 
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Figure 2: Our standard NNLO valence dsitributions together with their ilu uncertainties 
ai = 4 GeV^, according to the input parameters in Table 2 at Qq = 2 GeV^ for 
the central curves. Our standard NLO results j8] are shown by the dotted curves. For 
comparison the standard NNLO results of AMP06 [17] and BBG06 [87j are shown as well. 
Our dynamical valence distributions at = 4 GeV^ practically coincide with the standard 
ones shown. 
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Figure 3: Comparing the ±lo" error bands of our dynamical (dyn) and standard (std) 
NNLO gluon distributions at small x for various fixed values of Q^. Note that = 2 
GeV^ is the input scale of the standard fit. The central NLO results are taken from [8] 
with uncertainties comparable to the ones shown for NNLO for above 2 GeV^. For 
comparison the 'standard' NNLO results of AMP06 [T7] are shown as well. The results at 
= 2 and 20 GeV^ have been multiplied by 0.5 and 2, respectively, as indicated in the 
figure. 
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Figure 4: As in Fig. 3 but for the sea quark distribution x{u + d). 
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Figure 5: Comparison of our dynamical (dyn) and standard (std) NNLO small-x results 
for with HERA data for > 1.5 GeV^ [70-74]. The dynamical NLO results are 

taken from To ease the graphical presentation we have plotted F^^jQ^) + 0.5i{Q^) 
with i{Q^) indicated in parentheses in the figure for each fixed value of Q^. 
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Figure 7: Typical ilex uncertainty bands of our dynamical NNLO and NLO results in 
Fig. 5 for two representative values of . The 'standard' NNLO and NLO results are 
very similar. To ease the visibility we have added 0.1 to the results for = 8.5 GeV^ as 
indicated. 
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4 5 6 7 8 9 10 11 12 13 4 5 6 7 8 9 10 11 12 13 14 

M (GeV) 

Figure 8: Our dynamical NNLO and NLO [8j results, together with their ±la uncertainties, 
for Drell-Yan dilepton production in pp collisions for various selected average values of xf 
using the data sets of [HS]. For comparison the NNLO AMP06 results [.IT] are shown as 
well. To ease the graphical presentation we have multiplied the results for the cross sections 
by 10* with i indicated in parentheses in the figure for each fixed average value of xp- 
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Figure 9: As in Fig. 8 but for pd collisions. 
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Figure 10: Our dynamical NNLO and NLO results, together with their ±lcr uncertainties, 
for (j^'^/2(j^^ as a function of the average fractional momentum X2 of the target partons. 
The dynamical NLO results are taken from [8], and the NNLO AMP06 ones from [T7] . 
The data for the dimuon mass range 4.6 GeV < M < 12.9 GeV are from 
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Figure 11: The dynamical NNLO predictions, together with their ±lcr uncertainties, for 
the 'reduced' DIS cross section ar{x, Q"^) = F2 — {y"^ /Y^)Fl. The uncertainty bands of our 
previous dynamical NLO results (dashed curves) are very similar in size [8] as the ones 
shown for NNLO. The HERA data for some representative fixed values of are taken 
from [70-74]. 
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Figure 12: Our dynamical NNLO predictions for Orix^Q^^ but for different proton beam 
energies relevant for most recent HERA-Hl measurements [96]. The ilo" uncertainty 
bands are similar to the ones shown in Fig. 11. Notice that the curves terminate when 
y = \. 
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Figure 13: Dynamical parton model LO, NLO and NNLO predictions for Fl{x,Q'^) to- 
gether with the ±1(7 uncertainty bands at NLO and NNLO. The heavy charm (c) contri- 
butions at LO (short-dashed curves) and NLO (long-dashed curves) are shown as well. 
The results at = 100 GeV^ are multiplied by 0.5 as indicated. 
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Figure 15: Dynamical and standard NNLO and NLO predictions for Fl{x,Q'^) at a fixed 
value ofW = 276 GeV. The NLO(MS) results are taken from [8]. The (partly preliminary) 
HI data [721 ESI [Ml [IDS] are at fixed W ~ 276 GeV. The more recent HI data [96], which 
correspond to smaller values of W (larger x and Q^), are compatible with the indirectly 
determined data shown. 



43 



